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The isotopic-yield distributions and kinematic properties of fragments produced in transfer- 
induced fission of ^'^^Pu and fusion-induced fission of ^®°Cf, with 9 MeV and 45 MeV of excita¬ 
tion energy respectively, were measured in inverse kinematics with the spectrometer VAMOS. The 
kinematic properties of identified fission fragments allow to derive properties of the scission config¬ 
uration such as the distance between fragments, the total kinetic energy, the neutron multiplicity, 
the total excitation energy, and, for the first time, the proton- and neutron-number sharing during 
the emergence of the fragments. These properties of the scission point are studied as functions of 
the fragment atomic number. The correlation between these observables, gathered in one single 
experiment and for two different fissioning systems at different excitation energies, give valuable 
information for the understanding and modeling of the fission process. 
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I. INTRODUCTION 

Nuclear fission is a complex process where a num¬ 
ber of different properties of the nucleus influence the 
characteristics of the resulting splitting. In general, ex¬ 
perimental observations are limited to few observables, 
making it difficult to isolate the impact of specific nu¬ 
clear properties on the whole process. As an illustra¬ 
tion, the fission-fragment mass yields have been mea¬ 
sured intensively during the last decades, being in fact 
the result of the proton and neutron sharing between 
the two fragments. Isotopic yields have been accessible 
only for the light fragments P, 0| whereas techniques of 
delayed gamma-spectroscopy led to measures limited in 
range and in precision Q . Few years back, the measure of 
the complete distribution of the fragment atomic number 
was feasible when based on techniques involving inverse 
kinematics [3] , leading to new observations that were not 
expected from the study of mass distributions [13 . Re¬ 
cently, the use of inverse kinematics using transfer- and 
fusion-induced fission Q has given access to the isotopic 
distributions of fragments issued from well-defined fis- 
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sioning systems. In this scenario, the neutron excess (de¬ 
fined as the neutron-to-proton ratio) of each fragment 
over the whole distribution is a new observable that al¬ 
lows the investigation of the neutron and proton sharing 
between fragments during their emergence, the effect of 
shell structure at scission deformation, and the way both 
fragments share the total excitation energy. In order to 
reproduce the data, fission models need to make assump¬ 
tions on the scission configuration (deformation and nu¬ 
cleon sharing) and the partition of excitation energy be¬ 
tween both fragments. However, scission-point proper¬ 
ties are very difficult to derive experimentally because 
the fragment detection occurs after neutron evaporation, 
hiding the exact neutron numbers of the fragments. In 
general, the elongation at scission of the fissioning system 
is revealed by the measure of the total kinetic energy Q 
while the deformation of fragments and the sharing of 
the excitation energy is derived from neutron multiplic¬ 
ities measured either directly 1,113 or indirectly from 
the velocity and energy of both fragments measured in 
coincidence [HI- 

In the present paper, the kinematic properties of the 
fission fragments are used to trace back information on 
the original split at scission. Due to the advantages of in¬ 
verse kinematics, the average neutron excess of the frag¬ 
ments can be derived experimentally for the first time, in 
addition to the total kinetic energy and the neutron mul¬ 
tiplicities. These results put forward unique information 
on the resulting equilibrium of neutron and proton shar¬ 
ing in the elongation process up to the scission point. The 
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FIG. 1. Mean values of the fission fragment velocity spectra as a function of the neutron number of the elements produced in 
the fission of ^'*°Pu (left panels) and ^®°Cf (right panels). The error bars show the second momentum of the velocity spectra. 


distance between the fragments and their total excitation 
energy at scission are also determined and discussed. 


II. FISSION VELOCITIES 

The characteristics of the fragment distributions from 
fission of two different compound nuclei, ^'*°Pu and ^^°Cf, 
have been recently investigated Q. Both fissioning sys¬ 
tems were produced using a 130^111 impinging onto 
a target at 6.1 AMeV, an energy 10% above the 
Coulomb barrier. The ^'^°Pu fissioning system was pop¬ 
ulated in (^^C,^°Be) transfer reactions. The detection of 
the target-like nuclei allowed the selection of the differ¬ 
ent reaction channels and, with the measurement of the 
energy and angle of the target-like nuclei, the excitation 
energy produced in the transfer reaction is reconstructed 
event by event [l3. In this work, it is assumed that all 
the excitation energy produced in the reaction is carried 
by the beam-like fissioning system. This approximation 
has been validated in a recent investigation of the trans¬ 
fer process between 12 q Following this 

assumption, ^“^^Pu was produced with an excitation en¬ 
ergy distribution centered on E* ^ 9 MeV and a width of 
6 MeV. In the case of ^®°Cf, the fissioning system was the 
result of fusion reactions between the beam and the 
target. The detection of the target-like nuclei allowed 
the selection of the different reaction channels. The VA- 


MOS spectrometer mill was used for the identification 
of the fission fragments and permitted a measure of the 
complete isotopic distribution of the fission fragments for 
these well-defined fissioning systems. 

At the focal plane of the spectrometer, the trajectory 
of the transmitted fission fragments was determined us¬ 
ing two position measurements, giving access to the focal 
plane vertical and horizontal angles. The measurement 
of the time of flight and a trajectory-reconstruction pro¬ 
cedure, based on the diagonalization of the spectrome¬ 
ter transport matrix [l5l| . were used to reconstruct the 
fragment velocity vectors. From this knowledge and the 
information on the kinematics of the compound nucleus, 
it is possible to reconstruct the fission-fragment velocity 
vectors in the reference frame of the fissioning nucleus, 
which defines the fission velocity vector. In the present 
work, the fission velocity determination is improved with 
respect to the previous work the slowing down of 
the fission fragments in the target has been taken into 
account, whereas it was considered as negligible previ¬ 
ously. For each fragment, the velocity measured in the 
laboratory is then corrected for the energy loss follow¬ 
ing the prescription of Ref. [l^, in which the different 
parameters are adjusted by means of LISE-I--I- simula¬ 
tions [i3- In addition, the velocity distributions of each 
fragment have been corrected for transmission cuts (angle 
and ionic charge states) that modified the mean value of 
the velocity distribution by few percents. The resolution 
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FIG. 2. (Color online) Average fragment velocity (V) as a function of the frag ment atomic number Z produced in the fission of 
(a) and (b). Red lines follow the prescription of Wilkins et al. |19l] with constant deformation and neck parameters, 

and the masses deduced in Sec. IIIII (see text for details). The insets show the ratio between the measured data and the 
calculations. 


of the resulting fission velocities in the reference frame 
of the fissioning system depends on the resolution of the 
velocity and angle in the laboratory reference frame, the 
beam-energy straggling, and the angle and energy of the 
target-like product in the case of transfer-induced fission. 
Considering a resolution of 0.4% on the velo city measure¬ 
ment and an angular resolution of 5 mrad [l 8 | , the reso¬ 
lution on the resulting fission velocity is better than 2 %. 
At the limits of the VAMOS acceptance, heavier frag¬ 
ments are less transmitted and mass distributions suffer 
from cuts in their heavier part. These cuts also impact 
the velocity determination. The effect of the acceptance 
on the velocity, which can modify the value of the mea¬ 
surement in ~ 0.5%, is included in the error bars. The 
resulting velocity distributions in the reference frame of 
the fissioning nucleus C(A, Z) from fission of ^"^^Pu and 
^®°Cf are displayed in Fig. [H 

The average velocity {V) for each atomic number Z is 
calculated with the fragment yields Y{Z, A) as: 


{V){Z) 


EV(Z,A)-Y(Z,A) 

A 


Y.y{z,a) 

A 


( 1 ) 


They are displayed in Fig. [5] for both systems. The 
steady decrease with increasing Z, observed for both 
cases, reveals the momentum conservation that induces 
smaller velocities as the mass of the fragment increases. 
In general, the velocities measured in fission from ^^°Cf 
are between 2 % and 10 % larger than those of ^^°Pu due 
in part to the stronger Coulomb repulsion between the 
fragments, which are heavier in the case of ^^°Cf. The 
average velocity i^) is compared with the prescription 
of Wilkins et al. [l^ (see Eos. ITT] and IT^. with constant 
deformation and neck parameters. In this formulation, 
the mass associated to each fragment atomic number was 
from Eq. [51 as explained in the following section. In both 
fissioning systems, the deformation parameters were kept 


at a constant value of Pi = ^2 = 0.625 for the two frag¬ 
ments while the neck distance is set in order to match the 
values of (V) at symmetric splits. In the case of ^"*^°Pu, 
the neck is then d = 1.9 fm long, while for ^®°Cf it needs 
to be increased to d = 2.2 fm. The inset of Fig. IHa) 
shows clear deviations of the measured velocities from fis¬ 
sion of ^"'^°Pu with respect to those calculated with fixed 
deformation and neck parameters around Z ~ 52 and 
Z ^ 42. These are the signature of smaller deformation 
in the scission configuration appearing for specific pro¬ 
ton and neutron numbers, as will be discussed later. In 
the case of ^®°Cf, the deviations show that for asymmet¬ 
ric splits, the measured velocities are larger than those 
calculated with fixed Pi = P 2 = 0.625 and d = 2.2 fm. 


III. RECONSTRUCTION OF THE FRAGMENT 
MASS AT SCISSION 

The fission-fragment velocities contain information on 
important properties of the scission configuration such 
as deformation and masses of the nascent fragments. In 
the reference frame of the fissioning system and due to 
momentum conservation, the ratio of the two fragment 
velocities is equal to the inverse ratio of the initial masses 
M *2 on each fission event: 

Yi2i = Mi (2) 

U 272 Ml*’ ^ ^ 

with 7 = [1 — (U/c)^]“^/^. 

In the present experiment, the velocity of one frag¬ 
ment per fission event is measured. This is done over 
the complete fragment production, as shown in the pre¬ 
ceding section. In both investigated cases, the excitation 
energy of the fissioning systems is not enough to allow 
proton evaporation along the fission process, therefore, 
for a fissioning system of atomic number Zps , the average 
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FIG. 3. (Color online) The average neutron excess (N)IZ of fragments produced in the fission of (a) and ^®°Cf (b) as 

a function of their atomic number Z, is displayed before (black markers) and after (blue open markers) neutron evaporation. 
It is compared to the estimation of the GEF code before (black dashed line) and after (blue dashed-dotted line) evaporation. 
Red long-dashed lines correspond to the predictions of the LD-SP model. 


characteristics, including momentum, of fragments with 
atomic number Z\ are associated to the complementary 
atomic number Zi = Zps — -^i- Under these conditions, 
the momentum conservation can relate both fragments 
as: 

(Ui 7 iM*) = {U 272 M*), (3) 

where the observables with subscript I correspond to 
the fragments with atomic number Z\ and those with 
subscript 2 correspond to its complementary fragment 
Zi = Zps — Z\. The analysis of the measured fragments 
shows that (Vi^iM*) can be approximated to {Viji){M*) 
with a relative error smaller than 0.4%. The average val¬ 
ues of mass excess permit to express M* as the product 
of the number of nucleons and the unified atomic mass 
unit, M* uA*, with a deviation smaller than 0.01% 
in A*. In addition, a detailed simulation of the isotropic 
neutron evaporation by fragments, within ranges of ve¬ 
locity, mass, and Z as the ones measured, confirms that 
its effect on the fragment velocity is smaller than 0.15%. 
These approximations allows the use of the relation of 
Eq.[2]with the measured velocities averaged for each frag¬ 
ment: 

(U272) (411) ■ ^ ^ 

Assuming that no neutron evaporation occurs from sad¬ 
dle to scission nor from the neck at scission, and knowing 
the mass of the fissioning system Aps, the average masses 
at scission can be deduced: 

/4*\ — <Vi7l> 

(A*)= Aps-(At). 

In the case of ^^°Pu fissioning system, with an av¬ 
erage excitation energy of 9 MeV, no pre-scission neu¬ 
tron evaporation is considered; however, a contamina¬ 
tion of 20% from ^“^^Pu Q renders the final fissioning 


mass Aps = 240.2 for transfer-induced fission. In the 
case of ^®°Cf, the excitation energy is sufficient to pro¬ 
duce neutron evaporation before fission with a multiplic¬ 
ity that depends on the excitation energy and angular 
momentum induced in the reaction. In a similar re¬ 
action with beam on angular anisotropy has 

been used to determine the root-mean-squared angular 
momentum induced in the reaction at an energy close to 
the Coulomb barrier (2^. A value of 24Ii was determined 
for this reaction. Calculations based on the Bass mod- 
elisation [ 2 l|, adjusted to reproduce the angular reaction 
of i60-b238u data, give an estimation of 20h in the case 
of the present experiment. This angular momentum and 
an initial excitation energy of 45 MeV lead to an aver¬ 
age fissioning system of mass Aps = 249.6, conside ring 
the different fission probabilities obtained from GEF 
predictions, where first-chance fission happens in more 
than 68 % of the events. This result is confirmed by 
the measurement of pre-scission neutrons emitted from 
a slightly lighter compound nuclei with similar excita¬ 
tion energy and larger angular momentum, for which a 
multiplicity of 1 is found [ 2 ^. The uncertainty esti¬ 
mated on the fissioning-nucleus average mass is restricted 
to 0.3 mass units, considering possible small variations on 
angular momentum and the validity of the model. The 
pre-scission neutron evaporation also modifies the aver¬ 
age excitation energy of the fissioning system, estimated 
with the same code in E* ~ 42 MeV. 

The experimental neutron excess of the fission frag¬ 
ments is known to vary with the asymmetry of the fission 
split 0. In the framework of unchanged charge density, 
both fragments keep the same neutron excess as the fis¬ 
sioning nucleus. Experimentally, it has been observed 
that the light fragments are less neutron rich than this 
prediction, while the heavy fragments are more neutron 
rich. This is known as the charge polarization. The neu¬ 
tron excess at scission {N*)/Z is calculated with the av- 
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erage masses (A*) obtained in Eq. [5]as: 

(iV*)/Z(Z) = 

and it is shown in Fig. [3] with full symbols as a func¬ 
tion of the fragment Z for ^^°Pu and ^®°Cf. The de¬ 
duced {N*)/Z at scission is compared in this figure with 
the measured {N)/Z after neutron evaporation, shown in 
open symbols. In this case, the average (N) is obtained 
from the measured masses and fragment yields Q as: 

j:a-y{z,a) 

A 

The error bars in Fig. [3] include the statistical uncer¬ 
tainty, the estimated deviation due to acceptance cuts, 
the evaluated uncertainty of the fission velocity in the ref¬ 
erence frame of the fissioning system, and of the mass of 
the fissioning system. In the case of post-neutron (N) jZ 
for ^^°Pu, the large error bars around symmetry are pro¬ 
duced by the combined effect of low yields and contami¬ 
nation from other channels. It is also important to note 
the different nature of the calculation of both post- and 
pre-neutron neutron excess: the former is evaluated with 
the direct identification in {Z, A) of the fission fragments, 
while the latter is derived from a continuous observable 
such as the velocity. 

The resulting {N*)/Z for fission of ^^°Pu shows a step 
behaviour with a sudden increase that peaks around 
Z 50, paired with a minimum around Z ^ 44. The 
subsequent neutron evaporation does not attenuate the 
step structure, which is also observed in {N)jZ (open 
symbols in Fig. [3]). In the case of fission of ^®°Cf, {N*)/Z 
shows a smooth behaviour and a steady increase with 
Z, while neutron evaporation modifies significantly the 
trend. The fragments observed after neutron evaporation 
show a rather constant neutron excess, with no charge 
polarization. The smaller error bars observed in the frag¬ 
ments of ^®°Cf reflect the larger statistics in the popula¬ 
tion of this fissioning system compared to that of ^"'^°Pu, 
due to the cross sections of fusion-fission reactions be¬ 
ing about a factor 50 larger than that of transfer reac¬ 
tions (T^. [^. 

The deduced {N*)/Z can be compared with a simple 
calculation based on the scission-point model introduced 
by Wilkins et al. 0 , where the most probable combina¬ 
tion of neutron excess of both fragments is defined as the 
one that minimizes the total energy of the system. In the 
model used in the present work, only macroscopic ener¬ 
gies of the deformed fragments are considered, following 
the liquid-drop energy prescription of Myers and Swiate- 
cki [^, including Coulomb interaction and neglecting 
nuclear interaction between the two touching fragments. 
The results of this liquid-drop scission-point (LD-SP) cal¬ 
culation are displayed in Fig. [3] They show a steady in¬ 
crease of the neutron excess of scission fragments with 
their atomic number, reflecting the need to increase the 


number of neutrons for counterbalancing the increasing 
Coulomb repulsion inside the heavier fragments. This 
behavior is observed in experimental data from ^®°Cf, 
where the compound nucleus fissions at E* ^ 42 MeV 
of excitation energy. Still, the experimental neutron ex¬ 
cess deviates from this smooth behaviour with a sort of 
“S” shape that renders the heavier fragments slightly 
more neutron rich, and thus lighter fragments less neu¬ 
tron rich, than expected from equilibration of macro¬ 
scopic strengths. Stronger deviations are observed in the 
case of ^"^°Pu, which is produced with an average excita¬ 
tion energy of 9 MeV, due to the influence of the nuclear 
structure on the sharing of neutrons and protons between 
the fragments. Apart from these deviations, the deduced 
{N*)/Z and the calculation of the LD-SP model coincide 
for very asymmetric splits in both systems. 

A more realistic calculation, including the effect of nu¬ 
clear structure and neutron evaporation, can be done 
with the GEF code. Figure [3] shows the neutron excess 
calculated with GEF at scission and after neutron evap¬ 
oration for ^^°Pu and respectively. The structures 

appearing in the neutron excess of the ^^°Pu scission 
and final fragments are well reproduced by GEF, with 
a certain underestimation of the neutron excess for heav¬ 
ier fragments and therefore an overestimation for lighter 
ones. These structures are attributed to the effect of 
shells of protons and/or neutrons that tend to minimize 
the system energy around certain nucleon numbers in 
the nascent fragments, modifying their production when 
compared to pure macroscopic (or liquid-drop) behavior. 
At higher excitation energy, these effects are expected to 
be very much reduced (26[ and the calculations with the 
GEF code are consistent with this expectation. How¬ 
ever, a clear deviation of the experimental data with re¬ 
spect to the code appears in the case of ^^°Gf fission at 
E* ^ 42 MeV: while the calculations with GEF show a 
steady increase of the neutron excess with the Z, as ob¬ 
served in the data, a significant deviation of this slope 
is found between Z ~ 44 and Z ^ 54, rendering light 
fragments less neutron rich and heavy ones more neu¬ 
tron rich than those obtained with the code. Despite 
this difference, the GEF code reproduces the experimen¬ 
tal constant value of the measured {N)jZ of the final 
fragments as a function of their atomic number. The fea¬ 
tures of {N*)/Z and {N)jZ suggest the need to improve 
the description of the potential at scission, as well as the 
sharing of excitation energy, that modifies the neutron 
excess of fragments between scission and final fragments. 


IV. POST-SCISSION NEUTRON 
EVAPORATION 

The difference between the scission mass deduced from 
the velocity properties and the measured mass at the fo¬ 
cal plane of the spectrometer gives access to the average 
neutron multiplicity (z/) evaporated by the deformed scis- 
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FIG. 4. (Color online) Average neutron multiplicity (y) (a,b) and total neutron evaporation multiplicity (m+ 2 ) (c,d) for ^^°Pu 
and Blue dotted lines correspond to the calculations performed with the GEF code. 


sion fragments: 

{u){Z) = {A*){Z)-{A){Z). (8) 

The resulting average neutron multiplicities are shown 
in Fig. m In the case of ^'^°Pu, a saw-tooth behaviour of 
the neutron multiplicity, usually observed as a function 
of mass [i^, can be recognized despite the large fluctu¬ 
ations that appear in the tails of the function and for 
symmetric splits (see Sec. EH. The saw-tooth shape of 
the neutron multiplicity in low-energy fission is under¬ 
stood as the result of the deformation energy released 
when the scission fragments separate and recover their 
ground-state deformation 0- For higher-energy fission, 
the additional excitation energy is shared among the two 
fragments, while part of it might be evaporated during 
the saddle-to-scission path. This additional energy will 
be also released in part with evaporated neutrons, thus 
making the neutron multiplicity a probe of the energy 
sharing between the two fragments. In order to repro¬ 
duce the observed multiplicities, models describing the 
available phase space in each fragment on the level den¬ 
sity of a Fermi gas need to introduce additional parame¬ 
ters . A more accurate reproduction of the measured 
multiplicities is achieved with descriptions based on the 
evolution of the level density as the excitation energy 
increases from constant-temperature to Fermi-gas mod¬ 
els d^. The neutron multiplicities obtained with the 


GEF code, shown in Fig. E are calculated according this 
last description. The comparison with the data shows 
a similar behaviour, although with an overestimation of 
^ 1 neutron for light fragments and a similar underesti¬ 
mation for heavier fragments. These differences between 
the data and the code are related with those of {N*)/Z 
and (N) jZ previously described. The experimental mul¬ 
tiplicity iy) of fission fragments of ^^°Cf at E* ^ 42 MeV 
is displayed in Fig. Eb)- A steady increase with Z is 
observed, in agreement with direct measurements of neu¬ 
tron multiplicity in similar fusion-fission reactions 1,0. 
The kink observed at ^ ~ 54 can also be seen in data 
from (p-l-^'^^Pu) around A ^ 140 0. GEF calculations 
show a similar increase, though with a less steep slope. 
This difference may arise from the treatment on the evo¬ 
lution from constant-density to a Fermi-gas description 
and/or the description of the polarization at scission. 


It is worth noting that the comparison of the calcula¬ 
tions performed with GEF with the experimental data 
is more sensitive when studying the average multiplicity 
as a function of the fragment {y) than in the case of the 
post-scission total multiplicity per fission, defined as the 
sum of the multiplicities of fragment Z\ and its comple¬ 
mentary Zi = Z-ps — Zi, or the mass difference between 
the fissioning system and the fragments detected after 
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correspond to the prescription of Wilkins et al. (see Eq. El), displayed for reference, with constant deformation parameters 
Pi = P 2 = 0.625, and d — 1.9 fm for and d = 2.2 fm for ^®°Cf. Blue dashed lines correspond to the calculations of the 

GEF model. Asterisks show data from thermal-neutron induced fission of ^®®Pu [sj] . 


neutron evaporation: 

(*^1+2) = (j^i) + (1^2) = ^FS — (^1) — (^2)- ( 9 ) 

The total multiplicity ( 1 ^ 1 + 2 ) from fragments of ^^°Pu 
and g^]-g shown in Figs. |4l(c,d). The total multiplic¬ 
ity of ^'^°Pu shows an average value around (z^i-i- 2 ) ~ 3 
and an increase for symmetric and very asymmetric 
splits, where the deformation of fragments is expected 
to be maximal. However, as previously noted, the re¬ 
gion around symmetry for this system suffers from large 
uncertainties. In the case of ^^°Cf, a flat behaviour with¬ 
out structure is observed. The different features revealed 
in the study of {ly) are hidden in the case of the total 
multiplicity ( 1 / 1 + 2 )- In Fig. 21 {yi+ 2 ) is compared with 
calculations performed with the GEF code. As in the 
case of {v), the total multiplicity { 121 + 2 ) of ^"‘'’Pu is well 
reproduced, within the large statistical fluctuations. In 
the case of ^®°Cf, a similar constant behaviour around 
symmetry is observed in the GEF calculations. 


V. TOTAL KINETIC ENERGY AT SCISSION 

The characteristics of the total kinetic energy TKE 
distribution of fully accelerated fragments are mainly de¬ 
termined by the configuration of the system at scission 
and also, when measured in most experiments, by the 
following neutron evaporation. By using the fragment 
masses reconstructed at scission and their velocities, it 
is possible to derive the Coulomb potential energy that 
acts between the fragments that would be transformed 
in TKE once the fragments are fully accelerated. In or¬ 
der to simplify the nomenclature, this potential energy 
will be referred as average Total Kinetic Energy at scis¬ 
sion {TKE*), and it is the equivalent to the total kinetic 
energy after full acceleration but before prompt-neutron 
emission. This quantity keeps a similar information on 


the scission configuration as the TKE but without the 
influence of neutron evaporation. The average {TKE*) 
is calculated as: 

{TKE*) = u{Al) (( 71 ) - 1) + u{Al) (( 72 ) - 1), (10) 

with u = 931.494 MeV. The resulting {TKE*) is dis¬ 
played in Fig. [S] for both fissioning systems. Ignoring the 
contribution of pre-scission velocity, the total kinetic en¬ 
ergy at scission can be considered as the result of the 
Coulomb repulsion between the fragments at a distance 
D [H: 

TKE* = 1.44^32^. (11) 

the distance between the two charged centroids D can be 
expressed as [l^ : 

D = rQA\^^^ ^1 -I- EriiA^^^ ^1 -|- +d, (12) 

where ro = 1.16 fm. Pi and P 2 are the quadrupole de¬ 
formation parameters of fragments of mass number A\ 
and A2, and d is the distance between the surface of 
both fragments, associated with the length of the neck. 
Figure [5] shows a comparison between the experimental 
{TKE*) and this expression with constant parameters 
of deformation Pi = P2 = 0.625, and a neck parame¬ 
ter of d = 1.9 fm in the case of ^"'^°Pu and d = 2.2 fm 
in the case of ^®°Cf. In both cases, the neck parameter 
was chosen to match the experimental values at sym¬ 
metric splits, where the condition Pi = P 2 is trivially 
fulfilled and the fragment mass at scission is Afs/ 2. The 
masses of the fragments are those reconstructed at scis¬ 
sion, as explained in Sec. Ell The comparison with fis¬ 
sion from ^^°Pu shows strong deviations around Z 50 
and Z ^ 44. These values of {TKE*), higher by about 
10 MeV, correspond to more compact shapes, and are 
usually interpreted as the result of the superposition of 





different fission modes corresponding to the formation 
of fragments with atomic spherical and deformed closed 
shells Figure[5][a) shows the same features in thermal- 
neutron induced fission of ^^®Pu [3l| and a good agree¬ 
ment with calculations performed with the GEF code. 

In the case of ^^°Cf, the measured (TKE*) produces an 
almost flat behaviour around symmetric splits, between 
Z ~ 44 and Z ~ 54, while the difference with respect to 
the calculation with fixed deformation and neck parame¬ 
ters increases up to ~ 4 MeV towards asymmetric splits. 
The GEF code shows an overall underestimation from 
about ^ 4 MeV at symmetry to ^ 10 MeV for asymmet¬ 
ric splits, while it follows a similar shape as that with 
fixed scission parameters. In order to calculate {TKE*), 
GEF uses different approaches for fission with low- and 
high-excitation energy. In the former case, {TKE*) is 
calculated as the residue of the total energy after sub¬ 
tracting the collective, intrinsic, and deformation ener¬ 
gies. In the case of high excitation energy, GEF uses a 
similar approach as the LD-SP model, where the energy 
at scission is minimized to find the corresponding /3i and 
(32 deformation parameters, while keeping an empirical 
value of d = 3 fm for the neck. In addition, a certain 
amount of velocity pre-scission is added, and the influ¬ 
ence of fragment angular momenta is also considered. 

The relation between the elongation, namely, the dis¬ 
tance between fragments, and the resulting {TKE*) can 
be seen and quantified more clearly considering the ratio 
of the actual distance D between the fragments at scis¬ 
sion and their distance Dq as spherical touching nuclei. 
Following the prescription of Eq. [TTJ the average of this 
ratio {D/Dq) as a function of the fragment split can be 
calculated as: 


{D/Do) = 1-44 


Z1Z2 

{TKE*) 


jo{At)^^^+ro{A*)^ 


(13) 


Figure m shows {D/Do) for fission of both ^'^°Pu and 
250 cf The most compact shapes appear around Z 44 
and Z ~ 50, as observed in the {TKE*), while at sym¬ 
metric and very asymmetric splits, the elongation reaches 
values 60% larger than the distance between the corre¬ 
sponding spherical touching fragments. In the case of 
^®°Gf, the elongation is kept around a value of {D/Do) ^ 
1.60, with a slight increase at the symmetry and no ap¬ 
parent effects attributed to nuclear structure. Figures [5] 
and [6] show that the elongation of the system at scission 
changes with the asymmetry of the split, modifying the 
experimental {TKE*). Also in Fig. [51 the {D/Do) calcu¬ 
lated with the GEF code is shown. The comparison with 
the ^^°Pu data is consistent with the agreement found in 
both the neutron excess and the {TKE*). In the case 
qJ 250q£^ GEF predicts overall more elongated configura¬ 
tions with an almost constant value of D/Do ~ 1.66 and 
a slight dependence on Z somehow opposite to the one 
found in the experimental data: in the case of GEF, sym¬ 
metric splits are less elongated than asymmetric ones. 
Together with the situation found in the {TKE*), this 



FIG. 6. (Color online) Ratio between the distance between 
fragments at scission and that as spherical touching nuclei for 
fission of ^'*'^Pu (blue dots) and ^®°Cf (red squares). Dashed 
lines correspond to calculations done with the GEF code for 
^"'“Pu (blue dashed line) and ^®°Cf (red dashed line). 


behaviour suggests a more complex dependence on the 
fragment split of both the deformation and neck distance 
at scission with high excitation energies and relatively 
high angular momentum. 


VI. TOTAL EXCITATION ENERGY AT 
SCISSION 

The total energy available in the fission reaction E^ot 
is the sum of the mass of the fissioning system Mps in 
its ground state and its excitation energy Epg. At scis¬ 
sion, Etot is the sum of the masses of the fragments, 
the TKE*, and the excitation of collective (deformation, 
rotation, etc.) and intrinsic (single-particle excitations) 
degrees of freedom. The experimental determination of 
{TKE*) at scission allows then to determine the portion 
of the total energy that is transformed in the excitation 
energy of the fragments. Ignoring neutron evaporation 
during the saddle to scission path, the conservation of 
Etot links both instances: 

Etot = MfS + Apg . . 

= (M*) -b (M*) -b {TKE*) -b {TXE*), 

where {M*) is the average ground-state mass of the frag¬ 
ment i and {TXE*) is the total excitation energy avail¬ 
able for both fragments. The calculation of {M*) from 
the obtained {A*) is done by interpolation of the mass of 
the closest integer mass numbers: 

{M*)= {l-Wi)-mo{Z.,[{A*)\) + 

+Wi ■ Too {Zi, [(A*)J -b 1), (15) 

IV. = (a*)-L{A*)J, 

where toq {Zi, Ai) is the ground-state mass of the isotope 
{Zi,Ai) from the AME2012 evaluation [s^- The floor 
operator [(A*)J gives the largest integer not greater than 








9 





19 20 
0 


■a 

0 

N 


^ 10 - 
o 


.r i 


J_I_I_I_I_I_I_I_I_I_I_I_L 


1% I 


40 50 60 

Fragment Z 


FIG. 7. (a,b) Average difference between the ground states of the fissioning system and the fragments (Qgs) at scission (full 
symbols) and the l^KE*) (open symbols). (c,d) Total excitation energy at scission {TXEg^). (e,f) Post-neutron evaporation 
normalized yields for reference, (a), (c), and (e) correspond to while (b), (d), and (f) correspond to ^®°Cf. 


(A*). The quantity (TXE*) can also be regarded as 
the difference between the total energy Etot and the the 
potential surface at scission. By subtracting E^g, the 
depth of the potential can be expressed with respect to 
the ground state of the fissioning system as: 

(txe;,) =Mfs-{m^)-{m*)-{tke*) 

= (Qgs) - {TKE*), 

with (Qgs), displayed in Figs.[7Ka,b), as the average mass 
difference between the fissioning system and the frag¬ 
ments at scission. The error bars include the estimated 
deviation of the masses (M*) calculated from {A*) with 
Eq. [15] and those calculated from actual distributions of 
individual M*{Z,A). This deviation is found to be be¬ 
low 2 MeV. Both systems exhibit a (Qgs) with strong 
even-odd staggering, due to the systematic production 
of odd-odd or even-even Z pairs of fragments. The be¬ 
haviour of (Qgs) for each system is similar to that of their 
respective {TKE*) (see Fig.jS]), with two maxima around 
Z ^ 44 and Z ^ 50 in the case of ^'^°Pu and a flatter top 


in the case of ^^°Cf. Figures |7Kc,d) show the quantity 
(TXE*^) as a function of Z for ^^°Pu and ^®°Cf. Both 
cases have a similar behaviour with a well-defined maxi¬ 
mum centered at symmetric splits that raises from ~ 20 
up to ~ 25 MeV in ^^°Pu and from ~ 25 up to ^ 45 MeV 
in ^®°Cf. In Ref. [1^, a liquid-drop model (LDM) that 
considers surface, surface asymmetry, and Coulomb con¬ 
tributions is applied to thermal-neutron induced fission 
of 239p^, resulting in a TXE*^ = 15.27 MeV compared to 
~ 25 MeV determined in the present work. The same cal¬ 
culation applied to spontaneous fission of ^®^Cf results in 
TXE*^ = 24.94 MeV, much smaller than the ^ 45 MeV 
found for ^®°Cf. The main differences between the sys¬ 
tems used in the LDM calculations and those of this work 
correspond to the excitation energies and angular mo¬ 
menta produced in the reactions inducing fission. In the 
case of ^^'^Pu, thermal-neutron induced reactions yield 
an excitation energy of E* = 6.53 MeV and negligible 
angular momentum, while the transfer-induced fission of 
this work produces a distribution of excitation energy 
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FIG. 8. (Color online) (a) Minimum fraction of the total excitation energy {TXE*) evaporated with neutrons for ^'*°Pu (blue 
dots) and ^®°Cf (red squares), (b) Minimum energy released in neutron evaporation and excitation energy shared according 
with the mass ratio of the fragments for ^'*°Pu (blue dots and dashed line) and ^®°Cf (red squares and long-dashed line). 


centered at i?* ~ 9 MeV and a distribution of angular 
momentum of the order oi L ^ Qh (see [s^ for a simi¬ 
lar reaction). The spontaneous fission of ^®^Cf induces 
neither excitation energy or angular momentum, while 
the fusion-induced fission measured in this work produces 
^®°Cf with E* ^ 42 MeV and an average L ^ 20h. 

The presented values of (D/Dq) and {TXE*^) as func¬ 
tions of Z (Figs, m and [T]) can be interpreted as the shape 
of the potential landscape that the fissioning system ex¬ 
periences at scission. In both systems, symmetric splits 
are characterized by elongated shapes and maxima of po¬ 
tential. In the case of signatures of the most pro¬ 

duced asymmetric yields, around Z ^ 42 and Z ^ 52, 
are well correlated with shorter distances between the 
fragments and features in (TKE*) around Z ^ 44 and 
Z ^ 50 (see Fig.jSKa)), that can be interpreted as the ef¬ 
fect of nuclear structure. However, the deduced (TXE*^) 
for the same system does not show any distinct behaviour 
around these atomic numbers. The maxima in {TKE*) 
and (Qgs) compensate when both are subtracted to cal¬ 
culate (TXE*^), leaving no clear signature in the depth 
of the potential at scission. 

The total excitation energy gained at scission is re¬ 
leased by the fragments via post-scission neutron and 
gamma evaporation until their ground state is reached. 
The energy needed to evaporate these neutrons is the sum 
of their binding energies and their kinetic energy, which 
is a distribution between 0 and some MeV (see @,ll3 for 
instance). Therefore, the sum of the separation energies 
of the evaporated neutrons is a lower limit for the energy 
released by neutron evaporation. From the total neu¬ 
tron multiplicity {vi+ 2 ) (see Sec. IIVI) and the post- and 
pre-evaporation masses, this limit for the average energy 
released by neutron evaporation can be estimated 

as a function of the fragment split as: 

EZ \2 = {Mi) + (M 2 ) + m„ • { 121 + 2 ) - {M*) - (M*), (17) 

where (MJ and (M*) are the ground-state masses post- 
and pre-evaporation of fragment i, and m„ is the neutron 


mass. The neutron multiplicity evaporated after scission 
contains information on the excitation energy, including 
deformation energy, stored by the fragments at scission. 
In order to access this information from multiplicity mea¬ 
surements, is important to know the fraction of excita¬ 
tion energy carried by neutron evaporation. Figure [SKa) 
shows the ratio between i ?™(+2 and the deduced {TXE*) 
for the two systems studied. The behaviour of is 

similar to the ones observed in the neutron multiplicity 
in Fig. [4l as expected. In both systems, the fraction of 
{TXE*) evaporated with neutrons is smaller at symme¬ 
try, with ^^°Pu reaching some 50% and ^®°Cf showing a 
flat behaviour between Z ^ 44 and ^ 54 with more than 
60% of {TXE*). The two systems take more {TXE*) 
to evaporate neutrons, reaching up to 80%, as the split 
becomes more asymmetric, although a flatter behaviour 
would be within the error bars. 

The minimum energy released in neutron evaporation 
can be also calculated as a function of the fragment Z: 

E^^{Z) = {M){Z) + • {v){Z) - {M*){Z). (18) 

Figure [5Kb) shows as a function of Z for both sys¬ 
tems. As in the case of i?™i+ 2 , the behaviour is similar 
to the ones observed in the neutron multiplicity in Fig. 31 
with the heavy fragments, particularly in releasing 

more excitation energy in the form of neutron evapora¬ 
tion. This energy is a fraction of the excitation 

energy that each fragment stores up to scission, the rest 
being released as gamma emission. For fission with high 
E'pg, where structure effects vanish, it is expected that 
the total TXE is shared between the fragments accord¬ 
ing with their mass ratio at scission [s^: 

E:^+{Z) = {TXE*)^^^, (19) 

being E*^+ the excitation energy of the fragment, fol¬ 
lowing this prescription. It is important to note that, in 
any case, structure effects would not modify the E* pe¬ 
tition for symmetric or very asymmetric splits (see [36l |. 
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for instance). For lower energy fission, this behaviour is 
less valid, even at very asymmetric splits, for reasons dif¬ 
ferent from pure structure effects. Direct measurements 
of neutron multiplicity show that the heavy fragment 
keeps more excitation energy than Eq. [T9]or structure ef¬ 
fects predict [s^ . The evolution of the level density and 
the temperature with the mass of the fragment prevents 
them to reach equal temperatures, with the heavier frag¬ 
ment always cooler than the lighter one, and thus drag¬ 
ging more energy, resulting in an excitation energy of the 
hevier fragments larger than the predicted by Eq. mm. 
Figure [S^b) compares the deduced E™™ with In 

both cases, the E™™ of heavier fragments reaches similar 
values as their corresponding E*^^, while light fragments 
take less than the half of it. This behaviour suggests 
that, in both systems, the actual E* of the heavy frag¬ 
ments should be larger than E*^r^ in order to allow for the 
kinetic energy of the neutrons and gamma evaporation. 
This observation favors the description of a continuous 
flow of E* from the light to the heavy fragments due to 
their unbalanced temperatures. 


VII. CONCLUSIONS 

With the new generation of high-quality fission data 
obtained in inverse kinematics, new observables of the fis¬ 
sion process are available. From the reconstructed fission 
kinematics, the total kinetic energy (TKE*) and average 
masses of the two fragments can be reconstructed at scis¬ 
sion as a function of the fragment atomic number. In this 
work, the neutron excess of fragments is determined at 
scission for two different fissioning systems, ^"‘^Pu and 
^®°Cf, with different excitation energies. For the first 
time, an experimental insight on the proton and neutron 
sharing during the elongation process is put forward. In 
addition, the measurement of the total kinetic energy at 
scission allows the determination of the distance between 
fragments and the total excitation energy available. 

In the low-energy fission of ^^°Pu, {TKE*) and the 
mass difference (Qgs) show similar structures that may be 
a signature for the formation of little-deformed, closed- 


shell nuclei. These structures compensate when sub¬ 
tracted to estimate the total excitation energy (TXE*^) 
released in fission. This leads to an almost constant en¬ 
ergy release of about {TXE*^) ^25 MeV over the com¬ 
plete fragment distribution. The symmetric fission is ob¬ 
served to be associated with a slightly more energetic 
release of almost 30 MeV. In the higher-energy fusion- 
induced fission of ^®°Cf, with no neutron evaporation con¬ 
sidered between saddle and scission, the neutron excess 
of the scission fragments shows a steeper change around 
symmetry when compared to LD-SP and GEF models, 
resulting in more neutron-rich heavy fragments. The be¬ 
haviour of {TKE*) points to a smooth evolution of the 
scission configuration with the fragment split, with a flat 
behaviour around symmetry. This behaviour, together 
with the estimation of the distance between the frag¬ 
ments, indicates an almost constant elongation, though 
different from simple liquid-drop model considerations. 
Also in ^^°Cf higher-energy fission, a strong evolution of 
{TXE*g) is observed, with 25 MeV more excitation en¬ 
ergy in the symmetric splits than in the very asymmetric 
splits. 

A lower limit for the fraction oiTXE* released by neu¬ 
tron evaporation was found to evolve with the fragment 
split and to be above 50% of the total in both systems. 
The calculation of a lower limit of excitation energy re¬ 
leased by neutron evaporation as a function of Z sug¬ 
gests that the partition oiTXE* between the fragments 
according to their masses is not valid for these systems 
with Epg ~ 9 and ~ 42 MeV; being more suitable the de¬ 
scription with unbalanced temperatures and continuous 
flow of energy from the light to the heavy fragment. 
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